Two Particle States in an Asymmetric Box 



Xin Li a and Chuan Liu a 

^Department of Physics 
Peking University 

-^J- Beijing, 100871, P. R. China 

o ' 
o 



■ Abstract 

<D 

The exact two-particle energy eigenstates in an asymmetric rectangular box with 
periodic boundary conditions in all three directions are studied. Their relation with 
the elastic scattering phases of the two particles in the continuum are obtained. 
These results can be viewed as a generalization of the corresponding formulae in a 
cubic box obtained by Liischer before. In particular, the s-wave scattering length is 
(f-} \ related to the energy shift in the finite box. Possible applications of these formulae 

O ' are also discussed. 
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1 Introduction 



In a series of papers, Liischer obtained results [1,2,3,4] which relate the energy 
of a two-particle state in a cubic box (a torus) with the elastic scattering phases 
of the two particles in the continuum. This formula, now known as Liischer's 
formula, has been utilized in a number of applications, e.g. linear sigma model 
in the broken phase [5], and also in quenched QCD [6,7,8,9,10,11,12,13]. Due 
to limited numerical computational power, the s-wave scattering length, which 
is related to the scattering phase shift at vanishing relative three momentum, 
is mostly studied in hadron scattering using quenched approximation. CP- 
PACS collaboration calculated the scattering phases at non-zero momenta 
in pion-pion s-wave scattering in the 1 = 2 channel [12] using quenched 
Wilson fermions and recently also in two flavor full QCD [14]. In typical 
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lattice QCD calculations, if one would like to probe for physical informa- 
tion concerning two-particle states with non-zero relative three momentum, 
large lattices have to be used which usually requires enormous amount of 
computing resource. One of the reasons for this is the following. In a cu- 
bic box, the three momenta of a single particle are quantized according to: 
k = (27r/L)n = (27r/L)(ni, 712, 77,3), with n e Z 3 . 2 In order to control lattice 
artifacts due to these non-zero momentum modes, one needs to have large val- 
ues of L. One disadvantage of the cubic box is that the energy of a free particle 
with lowest non-zero momentum is degenerate. This means that the second 
lowest energy level of the particle with non- vanishing momentum corresponds 
to n = (1, 1, 0). If one would like to measure these states on the lattice, even 
larger values of L should be used. One way to remedy this is to use a three 
dimensional box whose shape is not cubic. If we use a generic rectangular box 
of size ijjiL) x {rj 2 L) x L with 7/x and 772 other than unity, we would have three 
different low-lying one-particle energy with non-zero momenta corresponding 
to n = (1, 0, 0), (0, 1, 0) and (0, 0, 1), respectively. This scenario is useful since 
it presents more available low momentum modes for a given lattice size, which 
is important in the study of hadron-hadron scattering phase shift. Similar sit- 
uation also occurs in the study of K to nir matrix element (see Ref. [15] for 
a review and references therein). There, one also needs to study two-particle 
states with non- vanishing relative three-momentum. Again, a cubic box yields 
too few available low-lying non-vanishing momenta and large value of L is 
needed to reach the physical interesting kinematic region. In these cases, one 
could try an asymmetric rectangular box with only one side being large while 
the other sides moderate. In an asymmetric rectangular box, the original for- 
mulae due to Luscher, which give the relation between the energy eigenvalues 
in the finite box and the continuum scattering phases, have to be modified 
accordingly. The purpose of this letter is to derive the equivalents of Liischer's 
formulae in the case of a generic rectangular (not necessarily cubic) box. 

We consider two-particle states in a box of size {i]iL) x (7/2-^) x L with periodic 
boundary conditions. For definiteness, we take 7/1 > 1,772 > 1, which amounts 
to denoting the the length of the smallest side of the rectangular box as L. 
The following derivation depends heavily on the previous results obtained in 
Ref. [3]. We will take over similar assumptions as in Ref. [3]. In particular, 
the relation between the energy eigenvalues and the scattering phases derived 
in the non-relativistic quantum mechanical model can be carried over to the 
case of relativistic, massive field theory under these assumptions, the same 
way as in the case of cubic box which was discussed in detail in Ref. [3]. For 
the quantum mechanical model, we assume that the range of the interaction, 
denoted by R, of the two-particle system is such that R < L/2. 

The modifications which have to be implemented, as compared with Ref. [3], 
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are mainly concerned with different symmetries of the box. In a cubic box, 
the representations of the rotational group are decomposed into irreducible 
representations of the cubic group. In a generic asymmetric box, the symmetry 
of the system is reduced. In the case of r\\ — rj 2 ^ 1, the basic group becomes 
D4; if 771 ^ 772 7^ 1, the symmetry is further reduced to D 2 , modulo parity 
operation. Therefore, the final expression relating the energy eigenvalues of 
the system and the scattering phases will be different. 



2 Energy eigenstates and singular periodic solutions of Helmholtz 
equation 



As discussed in Ref. [3], the energy eigenstates in a box is intimately related 
to the singular periodic solutions of the Helmholtz equation: 

(V 2 + k 2 )iP(r) = . (1) 



These solutions are periodic: ip(r + nL) = ^(r) and are bounded by certain 
powers of r = |r| near r = 0. The momentum modes in the rectangular 
box are quantized as: k = (27r/L)n. Here we introduce the notations: n = 
(n 1 /r) 1 ,n 2 /r]2,n 3 ) and n = (n^, r] 2 n 2 , n 3 ) with n = (ni,n 2 ,n 3 ) G Z 3 . For 
regular values of k, 3 the singular periodic solutions of Helmholtz equation 
can be obtained from the Green's function: 



We define: 3^m(r) = r l Yi m (fl r ), where fl r represents the solid angle parame- 
ters (9, 4>) of r in spherical coordinates; Y [m are the usual spherical harmonic 
functions. It is well-known that 3^ m (r) consist of all linear independent, ho- 
mogeneous functions in (x, y, z) of degree / that transform irreducibly under 
the rotational group. We then define: 

G lm (r;k 2 ) = y lm (V)G(v;k 2 ) , (3) 



which form a complete, linear independent set of functions of singular periodic 
solutions to Helmholtz equation. The functions Gi m (r; k 2 ) may be expanded 
into spherical harmonics: 



G lm (r;k 2 ) 



4:71 



Yi m (Q r )ni(kr) + J2 Mi m ;i'm'Y Vm >{Vt T )j v {kr) 



I'm' 



(4) 



This means that |fc| / (27r/L)|n| for any n € Z 3 . 
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Here, ji and ni are the usual spherical Bessel functions and the matrix M.i m ;i'm' 
is related to the modified zeta function via: 

M, m .„, = £ H,r% 1^; t ^V i^')(2fH)(2 ) + i) 





(z z' 




v o o oy 


m' 





(5) 



In this formula, the Wigner 3j-symbols can be related to the Clebcsh-Gordan 
coefficients in the usual way. For a given angular momentum cutoff A, the 
quantity M.i m ;Vm' can be viewed as the matrix element of a linear operator M 
in a vector space 7^a, which is spanned by all harmonic polynomials of degree 
I < A. The modified zeta function is formally defined by: 

Z lm (s,q 2 ; 771,772) = g (tf^ly ' ( 6 ) 

According to this definition, the modified zeta function at the right-hand side 
of Eq. (5) is formally divergent and needs to be analytically continued. Follow- 
ing similar discussions as in Ref. [3], one could obtain a finite expression for 
the modified zeta function which is suitable for numerical evaluation. It is also 
obvious from the symmetry of D4 or D 2 that, for I < 4, the only non-vanishing 
zeta functions at s = 1 are: Z 0Q , Z 20 , Z 2 ± 2 , Z^, Z^± 2 and Z±±±. One easily 
verifies that, if 771 = 772 = 1, all of the above definitions and formulae reduce 
to the those obtained in Ref. [3]. 

The energy eigenstates of the two-particle system may be expanded in terms 
of singular periodic solutions of Helmholtz equation. This solution in the re- 
gion where the interaction is vanishing can be expressed in terms of ordinary 
spherical Bessel functions, which is related to the scattering phases in the 
usual way. For the two-particle eigenstate in the symmetry sector T in a box 
of particular symmetry (either D 4 or D 2 ), the energy eigenvalue, E = k 2 /2fj, 
with n being the reduced mass of the two particles, is determined by: 

det[e 2l<5 - U(F)} = , U(T) = (M(T) + i)/(M{F) - 1) . (7) 

Here V denotes a particular representation of the group D 4 or D 2 . M(T) 
represents a linear operator in the vector space Ha(T). This vector space is 
spanned by all complex vectors whose components are t>/„, with I < A, and 
n runs from 1 to the number of occurrence of T in the decomposition of 
representation with angular momentum /, see Ref. [3] for details. To write out 
more explicit formulae, one therefore has to consider decompositions of the 
rotational group representations under appropriate symmetries. 
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3 Symmetry of an asymmetric box 



As mentioned in the beginning of this letter, modifications have to be made 
since the symmetry of an asymmetric box is different from that of a cubic 
one. We first describe the case rji = r\ 2 . The symmetry group is D 4 , which 
has 4 one-dimensional representations: Ai, A 2 , Bi, B 2 and a two-dimensional 
irreducible representation E. The representations of the rotational group are 
decomposed according to: 

= A+ , 1 = A 2 - + E-,2 = A+ + Bf + B+ + E + , (8) 



As is seen, in the Af sector, up to I < 2 both s-wave and rf-wave contribute. 
This corresponds to two linearly independent, homogeneous polynomials with 
degrees not more than 2, which are invariant under Z) 4 . These two polynomials 
can be identified as y 00 and y 2 o oc (x 2 + y 2 — 2z 2 ). Therefore, we can write 
out the reduced matrix element M.{Ai) n t = miy in this sector: 

^oo = VV 00 , m 20 = m 02 = - W 20 , 

m 22 = Woo + ^W 20 + j W 40 , (9) 
where we have introduced the notation: 

W Jm (l, q ; rji, V2) = tt 2 jti— • 10 



We find that, in the case of D± symmetry, Eq. (7) becomes: 

/ mL tan5 2 \ P 1 

1 + ^— tan5 = . 

\ m o 1 - m 22 tan d 2 J m 00 



Similar formula also appears in the case of cubic box except that the mixing 
with s-wave comes in at I — 4, not at I — 2. 

For the case rji ^ r/ 2 , the symmetry group becomes D 2 which has only 4 
one-dimensional representations: A, Bi, B 2 and B 3 . The decomposition (8) is 
replaced by: 

= A + , 1 = B^ + B 2 + £3- , 2 = A + + A + + B+ + B+ + B+ . (12) 

So, up to / < 2, A + occurs three times: once in I = and twice in / = 2. The 
corresponding basis polynomials can be taken as: 3^oo, 3^o and (3^2 +3^2-2 ) / V%- 
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If we denote the above three states as: 0, 2 and 2, the reduced matrix M(A + ) 
is three-dimensional with matrix elements moo, m 02 = 1^20 and 17122 given in 
Eq. (9) and the rest are given by: 



^02=^20 = --j=(W 2 2 + W2-2) , 



m 2 2 = rri22 



i^(W 22 + W2-2) + ^(^42 + W4-2) , 



2\/~E 1 / 5 

™22 = Woo - — W 20 + - W 40 + y ^( W 44 + W4-4) • 



(13) 



Similar to Eq. (11), the relation between the energy eigenvalue and the scat- 
tering phases now reads: 



m o 



cot 8q 



1 — m 2 2 tan 8 2 m 22 tan 8 2 
m 2 2 tan 8 2 1 — m^ tan 5 2 



(14) 



m 02 tan 2 8 2 
moo 



m 02 m 2 2 
m 2 (cot 8 2 - m22) 



m 2 tan 2 82 
m o 



m 2 (cot5 2 - m 22 ) 

^02 "^22 



If the ci-wave phase shift were small enough, it is easy to check that both 
Eq (11) and Eq (14) simplifies to: 



cot So(k) = moo — 



gooCjo g 2 ; yum) 

7i 3/2 r]ir] 2 q 



(15) 



For the general case, Eq. (11) and Eq. (14) offer the desired relation between 
the energy eigenvalues in the Af sector and the scattering phases for the cases 
771 = r] 2 and rji 7^ r] 2 , respectively. 



4 Large volume expansion of the scattering length 



It is known that in low-energy scattering processes, the scattering phases 8i(k) 
behaves like: tan<5;(/c) ~ k 2l+1 for small k, where k is the relative momentum 
of the two particles being scattered. It is easy to verify that both m o and 
m 02 behave like 1/q 3 as q ~ 0. Since tan<5 2 (g) goes to zero like q 5 , we see that 
the effects due to <i-wave phase shift in Eq. (11) and Eq. (14) are negligible as 
long as the relative momentum q is small enough. Therefore, in both cases, the 
s-wave scattering length ao will be determined by the zero momentum limit 
of Eq. (15). 
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For a large box, a large L expansion of the formulae can be deduced. Using 
Eq. (15), we find that the s-wave scattering length is related to the energy 
difference in a generic rectangular box via: 4 



SE = 



2ira 



Cl 



'■2 (//j, '/->) ( ^) + 



(16) 



Here, /x designates the reduced mass of the two particles whose mass values 
are mi and m 2 , respectively. Energy shift SE = E — mi — m 2 where E is the 
energy eigenvalue of the two-particle state. Functions 01(771,772) and 02(771,772) 
are given by: 



£00(1,05771,772) 
ci 771,772) = , 

7T77i772 

/ \ Z% (l,0;rii,r] 2 ) - Z 00 (2, 0; 771, 772) 

C2( " 1 '" 2)= ss 

where the subtracted zeta function is defined as: 



(17) 



(18) 



Table 1 

Numerical values for the subtracted zeta functions and the coefficients 01(771,772) 
and 02(7/1,7/2) under some typical topology. The three-dimensional rectangular box 
has a size Li = r/i-L, L 2 = r/ 2 L and L3 = L. 



Li : 


£2 


■L 3 


7/1 


% 


£00(1,05771,772) 


£ 00 (2, 0; 771,772) 


01(7/1,772) 


02(771,772) 


1 


: 1 : 


1 


1 


1 


-8.913633 


16.532316 


-2.837297 


6.375183 


6 


: 5 : 


4 


1.5 


1.25 


-12.964476 


41.526870 


-2.200918 


3.647224 


4 


: 3 : 


2 


2 


1.5 


-16.015122 


91.235227 


-1.699257 


1.860357 


3 


: 2 : 


2 


1.5 


1 


-10.974332 


32.259457 


-2.328826 


3.970732 


2 


: 1 : 


1 


2 


1 


-11.346631 


63.015304 


-1.805872 


1.664979 


3 


: 3 : 


2 


1.5 


1.5 


-14.430365 


53.784051 


-2.041479 


3.091200 


2 


: 2 . 


1 


2 


2 


-18.430516 


137.771800 


-1.466654 


1.278623 



4 For low relative momenta, the d-wave scattering phase behaves like: tan 5 2 (q) ~ 
a 2 k 5 = a 2 (2ir / X) 5 g 5 , with a 2 being the d-w&ve scattering length. If we treat the 
effects due to tan S 2 perturbatively, we see from Eq. (11) and Eq. (14) that, in 
Eq. (16), functions c\ and c 2 receive contributions that are proportional to (a 2 /L 5 ), 
which is of higher order in l/L for large L. 
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In Table 1, we have listed numerical values for the coefficients 01(771,772) and 
02(771,772) under some typical topology. In the first column of the table, we 
tabulated the ratio for the three sides of the box: 771 : 772 : 1. Note that for 
Vi — V2 — 1, these two functions reduce to the old numerical values for the 
cubic box which had been used in earlier scattering length calculations. 

5 Conclusions 

In this letter, we have studied two-particle scattering states in a generic rect- 
angular box with periodic boundary conditions. The relations of the energy 
eigenvalues and the scattering phases in the continuum are found. These can 
be viewed as a generalization of the well-known Liischer's formula. In particu- 
lar, we show that the s-wave scattering length is related to the energy shift by 
a simple formula, which is a direct generalization of the corresponding formula 
in the case of cubic box. We argued that this asymmetric topology might be 
useful in practice since it provides more available low-lying momentum modes 
in a finite box, which will be advantageous in the study of scattering phase 
shifts at non-zero three momenta in hadron-hadron scattering and possibly 
also in other applications. 
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